Spectral Components of the Radio Composite Supernova 

Remnant Gll.2-0.3 
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ABSTRACT 

We present a high-resolution radio study of the supernova remnant (SNR) Gil. 2— 0.3 
using archival VLA data. Spectral tomography is performed to determine the properties 
of this composite-type SNR's individual components, which have only recently been 
distinguished through X-ray observations. Our results indicate that the spectral index 
of the pulsar wind nebula (PWN), or plerion, is ap = 0.251q ?o- We observe a spectral 
index of as = 0.56 =L 0.02 throughout most of the SNR shell region, but also detect a 
gradient in a in the south-eastern component. We compare the spectral index and flux 
density with recent single-dish radio data of the source. Also, the radio efficiency and 
morphological properties of this PWN are found to be consistent with results for other 
known PWN systems. 

Subject headings: pulsars: individual (PSR J1811 — 1925) — supernovae: individual (Gil. 2— 0.3) 
— sunernova. remnants — radio continuum: TSM 



supernova remnants — radio continuum: ISM 

1. INTRODUCTION 

The proposed association between Gil. 2— 0.3 
and the supernova (SN) event of 386 AD 
(Clark & Stephenson 1977) has prompted 
many years of detailed study of this object. 
Radio and X-ray observations made with the 
Very Large Array ( VLA) and Einstein satel- 
lite by Downes (1984) verified the existence 
of a highly symmetric shell. Downes then 
argued that Gil. 2— 0.3 was likely the result 
of either a 300-500 yr old Type I SN explo- 
sion like Tycho or Kepler's supernova rem- 
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nant (SNR), or an explosion similar to the 
event which produced SNR Cassiopeia A, but 
which took place in 386 AD. Morsi &; Reich 
(1987) detected possible flat spectrum emis- 
sion in the central region from high-frequency 
radio maps made with the Effelsberg 100-m 
telescope, and concluded that Gil. 2— 0.3 be- 
longs to a plerion-shell composite SNR class. 
Overall spectral index estimates of a = 0.56 
(Downes 1984) and a = 0.49 (Morsi & Reich 
1987), where S v oc v a , were determined from 
old and recent integrated flux density mea- 
surements. Further VLA observations (Green 
et al. 1988) produced high resolution maps at 
20 and 6 cm, and revealed clumpy emission 
at the outer shell boundry, giving it the ap- 
pearance of an evolved Type II SNR similar 
to Cas A. Green et al. (1988) also estimated a 
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distance of ~ 5 kpc to the remnant based on 
its H I spectrum, and a diameter of ~ 6 pc. 
Reynolds et al. (1994) further supported the 
hypothesis linking Gil. 2— 0.3 to the historical 
event through correlations between radio and 
ROSAT X-ray observations, but argued for a 
Type la progenitor SN. 

Hard-spectrum, non-thermal X-ray emis- 
sion from a source within the remnant was 
detected in data from ASCA, which implied 
the existence of an embedded pulsar produc- 
ing plerionic emission (Vasisht et al. 1996), 
and confirmed that Gil. 2— 0.3 is a compos- 
ite remnant. A 65-ms X-ray pulsar (PSR 
J1811— 1925) was discovered by Torii et al. 
(1997); subsequent observations showed the 
characteristic spin-down age and rate of rota- 
tional kinetic energy loss to be r = 24, 000 yr 
and E = 6.4 x 10 36 erg/s, respectively (Torii 
et al. 1999). The exact location of pulsar PSR 
J1811— 1925 and the X-ray morphology of the 
pulsar wind nebula (PWN) were revealed for 
the first time in Chandra X-ray images (Kaspi 
et al. 2001). 

Recent high-frequency radio observations 
using the single-dish Effelsberg telescope re- 
vealed the spectral index of Gil. 2— 0.3 to be 
a = 0.50 for the combined spectrum of the 
PWN and SNR (Kothes & Reich 2001). A 
spectral index of a = 0.57 for the shell emis- 
sion alone was also calculated, by subtracting 
an estimated flat spectrum plerionic compo- 
nent. 

Because Gil. 2— 0.3 has long been mis- 
classified as a purely shell radio SNR, here 
we re-examine this source in light of its com- 
posite nature, as recently revealed by hard X- 
ray observations. In this paper, we re-analyse 
archival VLA data of Gil. 2— 0.3 to deter- 
mine its spectral index, using the technique of 
spectral tomography, which is a method that 
allows us to determine a for both the shell 
and plerionic regions separately. We com- 
pare the flux density and spectral parameters 



of each feature with results from single-dish 
radio studies. We also measure the structure 
and radio efficiency of the PWN to verify that 
they are consistant with other known plerionic 
remnants. These results will also be included 
in interpreting X-ray data results (Roberts et 
al., in preparation). 

2. OBSERVATIONS AND RESULTS 

Radio observations of Gil. 2— 0.3 were 
made with the VLA at 20 and 6 cm (L- and 
C-bands, respectively), between 1984 April 
and 1985 May. Details of these observations 
are summarized in Table 1. 

Data reduction and analysis were per- 
formed using standard procedures within the 
miriad package (Sault & Killeen 1999). The 
data were flux-density and antenna-gain cal- 
ibrated using the primary calibrator 3C286 
(J2000 1331+305). The secondary calibrators 
1749+096, 1748-253 and 1741-038 (B1950) 
were chosen for phase calibrations, based on 
observation coordinates and wavelength band. 
To aid in the process of detection and removal 
of bad visibilities, we performed calibrations 
on each data set individually. The header 
files were adjusted to match pointing centers; 
this step was necessary to prompt miriad 's 
restoring process into maintaining the small- 
est possible synthesized beam. 

Combined images at both 20 and 6 cm 
were formed using multi-frequency synthesis 
and a uniform weighting scheme. The images 
were deconvolved using a maximum entropy 
algorithm (Narayan & Nityananda 1986). We 
applied iterative phase self-calibration proce- 
dures in multi-frequency mode to make ad- 
ditional corrections to the antenna gains as 
a function of time. The resulting cleaned 
maps, as displayed in Figure 1, were convolved 
with a Gaussian restoring beam of dimensions 
276 x 178 and 27l x l'/5 for 20 and 6 cm, re- 
spectively. 
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Table 1 

vla observational parameters for gil. 2-0. 3 



Observing 


Array 


Frequencies 


Bandwidth 


Time on Source 


Date 


Configuration 


(MHz) 


(MHz) 


(min) 


1984 Apr 26 


u 


1465, 1515 


5U.UU 


39.5 


1984 Jun 06 


C 


1465, 1635 


50.00 


18 


1984 Jun 06 


C 


4535, 4965 


50.00 


19 


1984 Jul 05 


DnC 


1446, 1496 


12.50 


55.5 


1984 Jul 05 


DnC 


4923, 4873 


25.00 


144 


1984 Aug 27 


D 


4535, 4965 


50.00 


14.5 


1984 Sep 09 


D 


4835, 4885 


50.00 


34.5 


1985 Feb 22 


A 


1415, 1490 


6.25 


77 


1985 Feb 25 


A 


1415, 1490 


6.25 


82.2 


1985 May 19 


B 


1415, 1490 


6.25 


62.3 


1985 May 19 


B 


4640, 4670 


6.25 


129 


1985 May 25 


B 


1415, 1490 


6.25 


63.3 


1985 May 25 


B 


4640, 4670 


6.25 


112.8 



2.1. Spectral Index Determination 

To measure the spectral index a, where 
S u oc v~ a , from the 20- and 6-cm images, we 
first needed to spatially filter their u — v sky 
distributions in order to match spatial scales 
(Gaensler et al. 1999; Crawford et al. 2001). 
We degraded the nearly-complete u — v cov- 
erage of the 20-cm image to that at 6 cm 
through a process of u — v modelling, which 
essentially created a 20-cm visibility dataset 
with the same coverage as for the 6-cm image. 
This was imaged and deconvolved in the same 
way as were the original 6-cm data, and both 
maps were convolved to 10" resolution to aid 
clarity in the following tomography process. 

We could then compare the two images in 
order to make a spectral index determina- 
tion using the method of spectral tomography 
(Katz-Stone & Rudnick 1997). A difference 
image I at was calculated by scaling the 6-cm 
image by a trial spectral index at, and sub- 



tracting it from the 20-cm image, according 
to the formula 




where I20 and 1$ were the images being com- 
pared, and ^20 and v% were the average fre- 
quencies of the combined 20- and 6-cm im- 
ages, respectively. When the trial spectral 
index reached the actual spectral index of a 
particular feature, i.e. at = a, that region ap- 
peared to vanish into the local background of 
the difference image. Uncertainties in a were 
estimated by finding the range of at beyond 
which significant residuals appeared. If a cer- 
tain area of the feature possessed a slightly 
different spectral index, it appeared as a dis- 
tinctly positive or negative residual compared 
to the rest of the difference image. 
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Fig. 1. — High resolution images of SNR 
Gil. 2-0.3 at 20 (top) and 6 (bottom) cm. 
The cross indicates the location of its associ- 
ated pulsar PSR J1811-1925. Note that the 
size of the cross is much larger than the po- 
sitional uncertainty (Kaspi et al. 2001). The 
"stripes" in the 6-cm image are an artifact of 
emission from a bright point source located to 
the south (not shown). 



2.2. Tomography Results 

By examining the series of difference im- 
ages, we find approximate values of a for the 
PWN and SNR regions of Gil. 2-0.3. Posi- 
tive residuals appear as light emission against 
the neutral grey background, and negative as 
dark emission. 




Fig. 2. — Tomographic spectral index im- 
ages of the PWN region of Gil. 2-0. 3. The 
greyscale shows difference images for at = 
0.15 to 0.40 in steps of 0.05. Contours corre- 
spond to 6-cm data convolved with a 10" cir- 
cular beam (shown in lower right corners), at 
levels of 5% to 41% of peak flux 91 mJy/beam, 
in steps of 12%. 

Figure 2 shows a series of tomographic im- 
ages between 0.15 < at < 0.40 in Aa = 0.05 
intervals. The PWN region in the center dis- 
appears into the surrounding region within 
the shell when the trial spectral index is be- 
tween at ~ 0.20 and 0.30. Variations in the 
appearance of the shell interior between one 
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Fig. 3. — Tomographic spectral index images 
of the SNR shell region of Gil. 2-0. 3. The 
greyscale shows difference images for at = 
0.50 to 0.60 in steps of 0.02. Contours corre- 
spond to 6-cm data convolved with a 10" cir- 
cular beam (shown in lower right corners), at 
levels of 5% to 41% of peak flux 91 mJy/beam, 
in steps of 12%. 

image and the next suggest that the rms noise 
here roughly corresponds to a spectral in- 
dex difference of ~ 0.05, although we observe 
a very clear departure from uniform emis- 
sion between 0.15 and 0.30; consequently, we 
choose asymmetric uncertainties around our 
best estimate of ap = 0.25^q'io- 

In Figure 3, we show tomographic images 
between 0.50 < a t < 0.60 in Aa = 0.02 
intervals. The majority of the SNR shell 
fades into the surrounding background be- 
tween at ~ 0.54 and 0.58; we therefore es- 
timate a value of as = 0.56 ± 0.02 for most 
of the shell region. However, it appears that 
the SNR spectral index is not completely uni- 



form, most noticeably in the southeastern re- 
gion. At at = 0.52, dark patches can be seen 
in the outer shell; likewise, lingering positive 
residuals appear as lighter emission in the in- 
ner shell of the same region at at = 0.58. 
This gradient indicates that the spectrum for 
the outer shell is slightly flatter, and the in- 
ner shell slightly steeper, than for the bulk 
of the remnant. However, it should be noted 
that the similar size of the 6-cm beam relative 
to Gil. 2— 0.3 may have produced this appar- 
ent gradient if the primary beam is not per- 
fectly modeled. Mosaiced VLA observations 
are currently underway which will verify this 
feature. 

Our measured values for ap and as are 
in agreement with single-dish data as pub- 
lished by Kothes & Reich (2001), who predict 
a fitted spectral index of atotai = 

0.50 ± 0.02 

from integrated flux density values of the to- 
tal remnant. Assuming a flat spectrum core of 
ap = 0, they subtracted the pulsar wind com- 
ponent of the emission from the total, and in- 
ferred as = 0.57. This value is well within the 
uncertainty of our result. They also discuss 
shell structure, specifically, the southeastern 
region's appearance of being in a later evolu- 
tionary phase than the rest of the SNR, based 
on polarization properties and X-ray observa- 
tions. The possibility of a connection between 
their conclusion and our detected gradient in 
a should perhaps be investigated further. 

2.3. Flux Density 

Using MIRIAD, we determined the inte- 
grated flux over various regions in the source, 
and subtracted the background signal esti- 
mated from nearby sky. We estimate the 
background-subtracted flux density of the en- 
tire source to be St = 16.6 ± 0.9 Jy at 20 cm 
and 8.4 ± 0.9 Jy at 6 cm. Similarly, we sub- 
tract the background interior to the shell sur- 
rounding the plerionic region to find the PWN 
flux density S P = 0.36 ± 0.23 Jy (20 cm) and 
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0.32 ± 0.18 Jy (6 cm). The uncertainties are 
inferred from the rms variations of flux in the 
background; however, while the exterior back- 
ground was relatively uniform in brightness, 
the interior background was highly irregular, 
causing the large uncertainties in the PWN 
flux. We then estimate the SNR shell flux 
density to be S s = 16.2 ± 1.1 Jy (20 cm) and 
8.0 ± 1.1 Jy (6 cm) based on these measured 
results. 

The greatest source of systematic uncer- 
tainty in these values is a result of flux miss- 
ing due to the nature of interferometric mea- 
surements having incomplete u — v coverage. 
Single-dish measurements would not be sub- 
ject to the same uncertainty; consequently, in 
comparing our results with those of Kothes & 
Reich (2001), we estimate that ~ 5% of the 
total 20-cm flux and ~ 12% of 6-cm total flux 
is missing from our maps. We also measure 
the amount of missing 6-cm emission result- 
ing from incomplete coverage by comparing 
the flux levels found in the 20-cm image be- 
fore and after the u — v modelling process. 
This ~ 11% difference in the total object and 
~ 13% difference in the plerion indicates how 
much flux is lost in the 20-cm data by degrad- 
ing it to 6-cm visibility, thereby implying that 
the 6-cm data are also missing ~ 11%— 13% 
flux. 

3. DISCUSSION 

3.1. Spectral Properties 

We now compare the results of past studies 
(see Table 2) with our measurements of spec- 
tral index and flux density at 20 cm. To do 
this, we plot an integrated fitted spectrum of 
the form 

S T = S s (u/U6iy as +S P (u/U61)- ap , (1) 

where v is the frequency in MHz. Because of 
uncertainties in background subtraction and 
incomplete u — v coverage, we choose to plot 



an upper and a lower estimate for each value 
of flux density. We set S P = 0.36 ± 0.23 as 
our low value estimate, because we used a re- 
gion of high flux as background (i.e. inside 
the shell), and also because of missing inter- 
ferometric flux. By assuming the lowest pos- 
sible background (i.e. outside the shell), we 
find a high value for Sp, which gives a low 
limit estimate for Ss of 15.4 ± 0.9 Jy. In or- 
der to account for possible missing interfero- 
metric flux, we deduce high flux estimates by 
subtracting our lower bounded values from in- 
terpolated total power results of Kothes & Re- 
ich (2001). The high estimates for Ss and Sp 
are approximately 17.0 ± 0.9 Jy and 2.0 ± 1.6 
Jy, respectively. Figure 4 shows Equation 1 
plotted with our best high and low estimates 
of S, and maximum and minimum uncertain- 
ties of a. We find our values to be consistent 
with single-dish values, within estimated un- 
certainties. The higher frequency data points 
suggest that the true value of Sp may be close 
to our upper estimate; we therefore assume 
Sp ~ 1 Jy for our rough calculations of radio 
efficiency. 

3.2. Morphology and Age 

The tomography difference images allow 
us to identify each feature according to its 
spectral index. In particular, we recognize 
the region of central emission that belongs 
to the PWN as outlined in the as difference 
maps (Fig. 3) by the prominent area of dark 
emission. A contour of this region is pro- 
jected onto a detailed map of the plerionic 
region of Gil. 2— 0.3 in Figure 5. In this im- 
age, it is apparent that the protruding feature 
in the northwestern direction, and possibly a 
smaller, similar feature in the southeastern di- 
rection, is not spectrally part of the PWN, but 
may still be associated with emission from the 
interior of the remnant. A toroidal/jet mor- 
phology has been detected in the Crab PWN, 
so it is possible that we may also be seeing a 
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Table 2 

Integrated flux density measurements of Gil. 2-0. 3 



r rec[ utjncy 


H 111 v 1 lnnoifir 

r lux uuiihiLy 




(MHz) 


(Jy) 




330 


39.0 ± 2.0 


Kassim (1992) 


408 


34.8 ± 2.4 


Kassim (1989); Shaver & Goss (1970) 


1408 


17.7 ± 0.7 


Reich, Reich & Fuerst (1990) 


2695 


11.5 ± 0.5 


Reich et al. (1984) 


4850 


9.6 ± 0.5 


Kothes & Reich (2001) 


10450 


6.3 ± 0.4 


Kothes & Reich (2001) 


14700 


5.7 ± 0.4 


Kothes & Reich (2001) 


23000 


4.7 ± 0.5 


Kothes & Reich (2001) 


32000 


3.8 ± 0.4 


Kothes & Reich (2001) 




1000 10000 
Frequency [MHz] 



Fig. 4.— Log-log plot of Gil. 2-0.3 flux den- 
sity spectrum. See Table 2 for flux data 
points. The fits correspond to the power-law 
spectrum (Eq. 1), with the upper and lower 
bound parameters for Ss, Sp, as, and ctp, re- 
spectively, as listed in the legend. A detailed 
description of the upper and lower limits is 
found in § 3.1. 

similar shape here. If it is a torus, then it is 
most likely centered very near the pulsar po- 
sition, with an axis of ~ 30° west of north. 
There is some unexplained emission south of 
its eastern branch, which appears to be asso- 



ciated with the PWN but does not fit into the 
toroidal morphology, in addition to the pro- 
trusions to the NW and SE, which are unasso- 
ciated with the PWN. Further discussion on 
the plerionic structure will be presented by 
Roberts et al. (in preparation), as will the re- 
lationship between hard X-ray and radio emis- 
sion. 

We measure the angular diameter of the 
SNR to be ~ 4f5 at maximum extent, which 
corresponds to a linear diameter of D$ ~ 6.5 
pc, assuming the distance d ~ 5 kpc to the 
remnant is accurate (Green et al. 1988). To 
calculate the average diameter of the PWN, 
we find the total area encompassed by the ple- 
rion and approximate it as radially symmet- 
ric, resulting in an angular diameter of ~ lf3 
and Dp ~ 1.8 pc. The Chandra X-ray image 
shows a PWN diameter of nearly 1 f 3 (Kaspi 
et al. 2001), consistent with radio results. The 
ratio of diameters is ~ 0.28 in radio; we com- 
pare this value against results for other com- 
posite remants and find it to be higher than 
most other reported ratios (van der Swaluw & 
Wu 2001). Current theory predicts two main 
phases of PWN evolution inside the shell. In 
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3.3. Efficiency 




18 h ll m 34 5 32 3 30 s SB 5 26 3 

Right Ascension (J2000) 



Fig. 5. — Image of the PWN region in 
Gil. 2— 0.3 at 20 cm. The contour represents 
roughly the PWN region boundary, as indi- 
cated by the spectral tomography difference 
images. The cross indicates the location of 
the pulsar, as in Figure 1. 

the first phase, a bubble of plerionic gas freely 
expands into the supernova ejecta superson- 
ically, until it encounters the reverse shock 
of the SNR blast wave, which compresses the 
plerion. In the second phase, the PWN con- 
tinues to expand within the SNR at subsonic 
velocities, as described by a Sedov solution 
(Reynolds & Chevalier 1984; van der Swaluw 
& Wu 2001). Results of hydro dynamical nu- 
merical simulations predict the largest radial 
extent ratio at the point just before the re- 
verse shock has reached the PWN (Fig. 8 of 
van der Swaluw et al. 2000). We detect a large 
PWN radius relative to its SNR; therefore, 
this object has probably not yet encountered 
the reverse shock. Assuming these models 
are correct, this provides more evidence that 
Gil. 2— 0.3 is indeed a very young remnant, 
much younger than is indicated by the 24,000 
yr characteristic age of the pulsar (Torii et al. 
1999; Kaspi et al. 2001). 



Gaensler et al. (2000) characterize the effi- 
ciency e of PWN radio emission by the frac- 
tional contribution of spin-down energy to ra- 
dio luminosity, 

Lr 
e = — , 
E 

where Lr is the radio luminosity of the PWN 
and E is the spin-down luminosity of the asso- 
ciated pulsar. They establish upper limits on 
e by re-examining results predicted by Frail 
& Scharringhausen (1997), with more up-to- 
date estimates for interstellar medium densi- 
ties. Six of the eight other known radio PWN 
are reported to have a "typical" efficiency of 
e ~ 10~ 4 . Using our result of ap « 0.25 and 
integrating over the range of 10 MHz to 100 
GHz, we find 

VI Jyy \5 kpcj 

Thus, the efficiency of Gil. 2— 0.3 is consistent 
with those of other detected PWN. 

4. CONCLUSIONS 

We have imaged SNR Gil. 2-0.3 at high 
resolution and revealed the spectral properties 
of its individual components. Our spectral 
index and flux density results are consistent 
with recent studies of this source and other 
composite-type SNR. However, measurements 
of its morphology suggest that it is in an ear- 
lier phase of evolution (just entering the Se- 
dov phase) than most other detected PWN 
inside remnant shells, according to numerical 
models. This result provides more evidence 
linking Gil. 2— 0.3 to the supernova event of 
386 AD. We also find the radio efficiency to 
be as expected for an object surrounding a 
young, energetic pulsar. 

We wish to acknowledge Bryan Gaensler 
for helpful discussions and David Green for 
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